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ABSTRACT 

We have discovered an x-ray source, SAX J0635+0533, with a hard spectrum within 
the error box of the GeV gamma-ray source in Monoceros, 2EG J0635+0521. The 
unabsorbed flux from the source is 1.2 x 10~^^ ergcm"^ s~^ in the 2-10 keV band. The 
x-ray spectrum is consistent with a simple powerlaw model with absorption. The photon 
index is 1.50 ± 0.08 and we detect emission out to 40 keV. Optical observations identify 
a counterpart with a V-magnitude of 12.8. The counterpart has broad emission lines 
and the colors of an early B type star. If the identification of the x-ray/optical source 
with the gamma-ray source is correct, then the source would be a gamma-ray emitting 
x-ray binary. 

Subject headings: pulsars: general — gamma rays: observations — stars: individual 
(2EG J0635+0521, SAX J0635+0533, LSI 61° 303, 2CG 135+01, PSR J1259-63) — 
stars: neutron 
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1. Introduction 
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The unidentified gamma-ray point sources concen- 
trated along tlic Galactic plane present a long out- 
standing puzzle (Swanenburg et al. 1981). Their po- 
sitions, timing, and energy spectra should provide a 
means for their identification. However, the gamma- 
ray error boxes are sufficiently large that multiple op- 
tical and x-ray candidates are usually found for Galac- 
tic sources. The only convincing identifications have 
been made through detection of a pulsed gamma-ray 
signal at a period known from radio or x-ray obser- 
vations (Bignami & Hermsen 1983, Halpern & Holt 
1992, Thompson et al. 1994, Ramanamurthy et al. 
1995). The known gamma-ray pulsars have spectra 
significantly harder than those of most of the uniden- 
tified sources. Recently, Merck et al. (1996) stud- 
ied the energy spectra of the unidentified Galactic 
gamma-ray sources and identified only eight sources 
with energy spectra sufficiently hard to be pulsar can- 
didates. 

Here, wc report on x-ray and optical observa- 
tions in the field of one of the unidentified Galactic 
gamma-ray sources with a hard gamma-ray spectrum, 
2EG J0635-I-0521. In §2, wc report the discovery of an 
x-ray point source with an unusually hard spectrum 
within the error box of 2EG J0635-h0521. In §§3 and 
4, we describe the properties of an optical counter- 
part of the x-ray source. We conclude, in §5, with 
comments on identification of the x-ray and gamma- 
ray sources and on the possible nature of the source. 

2. X-ray observations 

Using the Satellite Italiano per Astronomia X (Bep- 
poSAX), we observed a field centered on a region of 
unusually high x-ray hardness and brightness found in 
an x-ray survey of the Monoceros supernova remnant 
(Leahy, Naranan, & Singh 1986) and located within 
the error box of the EGRET source 2EG J0635-F0521 
(Thompson et al. 1995). The observation was car- 
ried out from 23 Oct 1997 3:20:26 UT to 24 Oct 1997 
0:11:52 UT. The field of view of the BeppoSAX nar- 
row field instruments covered approximately half of 
the error box of 2EG J0635+0521. The x-ray im- 
ages from the LEGS and the MEGS 2 and 3 detec- 
tors reveal a single point source located at RA = 
06*^35'"17".4, DEG = +05°33'20".9 (J2000). There is 
a 1' systematic uncertainty in the position reconstruc- 
tion. We denote the source as SAX J0635-F0533. This 
source lies within the error box of 2EG J0635-I-0521 
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Fig. 1.— X-ray image of SAX J0635-I-0533 for the 

5 12 keV band. Shown are the position of the opti- 
cal counterpart (cross), and the 95% confidence error 
elHpse for 2EJ 0635-1-0521 (dashed line). 

(see Figure 1) and near the peak of the gamma-ray 
emission as reconstructed by Jaffe et al. (1997). The 
image is consistent with that of a single point source, 
see Fig. 1. We find no evidence for diffuse hard x-ray 
emission. 

The x-ray spectrum of the source is quite hard with 
a spectral index near 1.5 and significant emission de- 
tected out to 40 keV. Shown in Figure 2 are data from 
the LEGS, the combined MEGS 2 and 3, and the PDS. 
We used an extraction radius of 2' for the MEGS and 
4' for the LEGS, the standard background files, and 
response matrixes calculated for the source position 3' 
off-axis. The x-ray spectrum from all the instruments 
can be fitted with a simple powerlaw model with ab- 
sorption. This fit is acceptable with x^/DoF = 97/81. 
The fitted column density is (2.0 ± 0.3) x 10^^ cm'^ 
and the powerlaw index is 1.50 ± 0.08. The unab- 
sorbed fiux from the source is 1.2 x 10^^^ erg cm~^ s~^ 
in the 2-10 keV band. We also fitted an exponen- 
tially cutoff powerlaw which gave x^/DoF = 96/80. 
The exponential cutoff does not improve the fit. The 
90% confidence lower bound on the cutoff energy is 
37 keV. We also fitted a Raymond-Smith model spec- 
trum which gave a worse fit x^/DoF = 122/83 and 
an unphysically high temperature, kT > 50 keV. 

We searched for pulsations in the x-ray data over 
a period range from 0.030 s to 1000 s. We found 
no significant pulsed signal. From the period with 
the highest single trial significance, we place an upper 
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Fig. 2. X-ray spectrum of SAX J0635+0533. 
Shown are data from the LEGS (1.0-4.0 keV), MEGS 
(1.8-10 keV), and PDS (12-60 keV). The soUd line is 
a fit of a power law of photon index of 1.5 with ab- 
sorption. 

limit on the pulsed fraction of 26%. The continuum 
of the timing power spectrum is consistent with the 
Poisson level for frequencies above 0.002 Hz. 

This field was previously viewed with the Einstein 
IPC (region 3 in Leahy et al. 1986). The source was 
not detected with Einstein because, we now know, 
it fell on a strongback arm. This field was also ob- 
served for 4371 s by the ROSAT PSPC on 1992 Octo- 
ber 9. The position of SAX J0635-h0533 falls 34' off 
axis in this image, where the point-spread function 
is significantly degraded. Nevertheless, there may be 
a very weak source present at this position. Fold- 
ing the spectrum of SAX J0635-I-0533 through the 
ROSAT response, we find that the predicted number 
of counts is roughly consistent with the number ob- 
served. There is no strong evidence for variability of 
the source. 

3. Optical observations 

In order to identify the optical counterpart of 
SAX J0635-I-0533, we obtained optical spectra of 
seven stars within and around the 1' radius error cir- 
cle on 1997 Dec. 10, using the 2.4 m Hiltner tele- 
scope of the MDM Observatory. The spectra covered 
the wavelength range 4500-7500 A at 4 A resolution. 




Fig. 3.— Finding chart for SAX J0635-F0533 from 
the Digitized Sky Survey. Field is 4' x 4'. North is 
up, cast is loft. The X-ray error circle has a radius of 
1'. The position of the optical counterpart is (J2000) 
RA = 06^35^18^.29, Dec = -F05°33'6".3. 

The brightest star within the error circle (Figure 3) is 
the only one with broad emission lines and other spec- 
tral features that are characteristic of a reddened Be 
star. Be stars are often hard X-ray sources because 
of emission from an unseen compact companion such 
as a neutron star. No other star in the error box has 
optical spectroscopic properties resembling those of 
known hard X-ray sources. On this basis we identify 
the Be star with SAX J0635-I-0533. This star is listed 
in the USNO Al.O astrometric catalog of the Palomar 
Observatory Sky Survey (Monet et al. 1996), at posi- 
tion (J2000) RA = 06'>35'"18«.29, Dec = -F05°33'6".3. 
The B1950 coordinates are RA = 06^32'"38".ll, Dec 
= -|-05°35'34".3. The same catalog gives approximate 
magnitudes of B = 12.7 and R = 12.7. There is 
no 20 cm radio source present at this position in the 
NVSS catalog (Gondon et al. 1998). 

Additional optical spectra were obtained with the 
Goldcam spectrograph on the KPNO 2.1 m telescope 
on 1998 January 28, 29, and 30. On the first two 
nights the spectra covered the wavelength range 3867- 
7516 A with a resolution of 5-6 A, while on the third 
night the spectra covered the range 5548 -8613 A with 
a resolution of 4 A. Three sequences of 300 s exposures 
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Fig. 4. — Composite optical spectrum of SAX J0635+0533 obtained on the KPNO 2.1m telescope on 1998 January 
29 and 30. In the middle panel, the flux scale refers to the lower trace. The upper trace is the same spectrum 
multiplied by a factor of 6. 



were taken on each night (at the beginning middle, 
and end of the night), yielding a total of 39 spectra. 
Lastly, we obtained one more spectrum on the Hiltner 
telescope on 1998 April 7. 

All of the spectra obtained are consistent with each 
other in continuum flux and emission-line properties. 
The mean of several of the KPNO spectra is shown in 
Figure 4, and their emission-line properties are listed 
in Table 1. The quoted line widths have been cor- 
rected for instrumental resolution. The strong Ha 
emission line has a narrow core and very broad wings. 
It can be fitted well by the sum of a broad and a nar- 
row Gaussian. The narrow component has an equiva- 
lent width (EW) of 27 A and a full-width at half max- 
imum (FWHM) of 330 km s~^, while the broad com- 
ponent has EW = 8 A and FWHM = 1300 km s'^ 
The wings of Ha extend to a full width at zero inten- 
sity (FWZI) of 2100 ± 200 km s"!. 

A 300 s exposure on SAX J0635-I-0533 was ob- 
tained by Goeff Marcy using the HIRES instrument 
(Vogt et al. 1994) on the Keck I telescope, on 1998 
January 26. The resolution is 2 km s^^. Figure 5 
shows the profiles of H/3 and He I A5876, after normal- 
izing them by the local continuum level and smooth- 
ing with a 14 km running boxcar filter. The 



emission lines are double peaked, with peak sepa- 
ration of « 200 km s-i and FWZI « 500 km s'^ 
Thus it appears that the narrow component identified 
in the moderate-resolution spectra as having FWHM 
« 330 km s^^ is actually resolved into a disk- like pro- 
file, similar to what is seen in the Be star X-ray tran- 
sient A0535-F26 by Clark et al. (1998), and attributed 
by them to a circumstellar disk around the Be star. 

The individual KPNO spectra were used to search 
for radial velocity variations of the Ha line over the 
three nights, but no significant variations were found. 
The scatter in the measured velocities appears to be 
random, and defines an upper limit to the real radial 
velocity variations on this time scale of 30 km s~^. 

4. Mean Spectrum, Reddening, and Distance 

The observed Johnson magnitudes (Johnson & Mor- 
gan 1951) measured from the mean spectrum of Fig- 
ure 4 are B = 13.81, V = 12.83, R = 11.98. Ab- 
sorption in many diffuse interstellar bands (DIBs) is 
apparent in Figure 4, and the equivalent widths of 
these can be used to estimate reddening and distance. 
We measured four of the best-studied interstellar ab- 
sorption features; their EWs are listed in Table 1. In 
addition to the EWs, the depth of the A4430 feature 
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Fig. 5. — Keck HIRES spectrum of selected emission 
lines, H/3 and He I A5876, normalized to the contin- 
uum and smoothed with a 14 km running boxcar 
filter. 

was also measured to be 13% below the continuum 
(this quantity is a reddening indicator in addition to 
the EW). There are several calibrations of the DIBs 
in the literature. The most extensive work is that of 
Herbig (1975) who studied the correlation of the EWs 
of many DIBs with reddening. Wamplcr (1963, 1966) 
noted that the calibration of the A4430 feature as a 
reddening indicator depends on Galactic longitude. 
More recently, Tiig & Schmidt-Kaler (1981) also in- 
vestigated this issue and came up with their own cal- 
ibration for the depth of the A4430 feature. We esti- 
mated the extinction E{B -V) to SAX J0635-F0533 
using each of these calibrations, and adopted a mean 
value of E{B — V) = 1.2 ± 0.2, where the uncertainty 
corresponds to the dispersion among different cali- 
brations. Assuming that Ay = 3.2 E{B — V), we find 
Ay = 3.8 ± 0.6. 

The Mon 0B2 association is a cluster of massive 
stars at a distance of 1620 pc (Garmany & Stencel 
1992), and within 1° of SAX J0635-F0533. The H I 
column densities to 11 stars in this cluster have been 
measured with lUE using Lyman a absorption and 
are in the range 1.4-3.8 x 10^^ cm~^ (Diplas & Savage 
1984). The estimated E{B - V) to SAX J0635-F0533 
corresponds to Nu = 6.0 x 10^^ cm~^ using standard 
conversions (Savage & Mathis 1979). Thus we can 
be confident that SAX J0635-I-0533 is no closer than 
the 1620 pc distance to Mon 0B2. The mean column 



density to the 11 stars is 1.93 x 10^"'^ cm ^. Assuming 
that = 1 for N{H I) = 1.6 x lO^^cm'^, this 
corresponds to Ay = 0.8 mag/kpc. If we assume 
that 0.8 mag/kpc is a typical value along the line of 
sight, the distance to SAX J0635-I-533 is probably not 
further than 5 kpc. 

The total Galactic 21 cm column density in this 
direction is 7.0 x lO^^ cm'^ (Stark et al. 1992), 
which is larger than the column density estimated for 
SAX J0635-f-0533 from interstellar absorption. Thus, 
we regard the estimated extinction from the DIBs as 
compatible with the 21 cm H I column, while the 
larger X-ray measured value, TVh = (2.0 ± 0.3) x 
10^^ cm~^, implies that there must be some circum- 
stellar gas that is responsible for the extra X-ray 
absorption. We note that the omission-lino ratio of 
Ha/H/3 can be used with recombination theory to de- 
rive an upper limit to E{B — V) of 1.4, and that this 
upper limit is also well below the X-ray measured col- 
umn. 

After correction for reddening, the Johnson magni- 
tudes estimated from the spectrum become B = 8.74, 
V = 8.99, and R = 9.14. Such colors are character- 
istic of an early B type star. For a distance range of 
2.5-5 kpc, the absolute visual magnitude would be in 
the range —3.1 to --4.5, and the stellar classification 
in the range B2V to BlIII. 

5. Discussion 

Assuming 200 Be/X-ray binaries within 5 kpc 
(Rappaport & van den Heuvel 1982) distributed uni- 
formly in Galactic longitude and in a latitude dis- 
tribution with an rms width of 2.2°, wo estimate 
there is a 4% chance of finding a spurious overlap 
of a Be/X-ray binary in an error box the size of that 
for 2EG J0635-F0521. This probabiUty is reasonably 
low. In addition, both Be/X-ray binaries and the low- 
latitude EGRET sources are preferentially found near 
OB associations (Kaaret & Cottam 1996). This may 
suggest that some of the unidentified EGRET sources 
are, in fact, associated with Be/X-ray binaries. 

While identification of SAX J0635-h0533 and its 
optical counterpart with 2EG J0635-I-0521 must re- 
main tentative awaiting an improved gamma-ray po- 
sition or detection of periodicity in both x-rays and 
gamma-rays, the positional coincidence and the hard 
spectrum of the x-ray source are suggestive and here 
we entertain some speculations based on the possi- 
ble identification. Assuming the association of the 
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x-ray source SAX J0635+0533 with the optical coun- 
terpart discussed above and with 2EG J0635-I-0521, 
allows us to find the luminosity of the source based 
on the optical distance estimate. The 2-10 keV X- 
ray luminosity would be 9 — 35 x 10'^^ ergs s~^ and 
the gamma-ray luminosity for 0.1-5 GeV would be 
1.3 — 5 X 10^^ ergss^^, both assuming isotropic emis- 
sion. The x-ray to gamma-ray flux ratio is 0.07. If 
the identification is correct, then the source would be 
a gamma-ray emitting x-ray binary. 

It is possible that the binary contains a neutron 
star and that both the x-ray and gamma-ray emis- 
sion are magnetospheric pulsar emission. If the bi- 
nary is wide, as is tentatively suggested by our non- 
detection of radial velocity variations, then the pul- 
sar emission may be be relatively uninfluenced by the 
presence of the companion star. There are several 
pieces of evidence in favor of pulsar emission. The ra- 
tio of x-ray to gamma-ray flux is consistent with those 
of "Vela-like" pulsars. The power-law x-ray spectrum 
with photon index near 1.5 is similar to that of most 
isolated gamma-ray pulsars above 2 keV (Wang et 
al. 1998). The gamma-ray flux is constant over 14 
different observations (McLaughlin et al. 1996), as 
would be expected for a pulsar. And finally, as noted 
above, 2EG J0635-h0521 was identified by Merck et 
al. (1996) as one of only eight unidentified Galactic 
gamma-ray sources with energy spectra sufficiently 
hard to be pulsar candidates. 

Alternatively, either the x-ray emission, gamma- 
ray emission, or both may arise from interaction of the 
energetic particle wind from the pulsar with the wind 
or the radiation from the companion. Such mecha- 
nisms have been suggested for Be/millisecond radio 
pulsar binary PSR J1259-63 and the variable radio 
Be/X-ray binary system LSI 61° 303. In PSR J1259- 
63, emission up to 200 keV is thought to arise from a 
shock interaction of the energetic particle wind from 
a pulsar with the wind from the Be-star (Grove et al. 
1995). The x-ray photon index, 1.5-1.9 (Kaspi et al. 
1995), is similar to that we find for SAX J0635+0533. 
However, PSR J1259-63 was not detected by EGRET 
(Tavani et al. 1996) and the x-ray emission is highly 
variable. It is possible that the shock acceleration 
mechanism could operate to higher energies if the 
shock region is less compact than that of PSR J1259- 
63 at periastron (Tavani & Arons 1997). 

The gamma-ray source 2CG 135-1-01 has long had a 
tentative association with the unusually variable ra- 
dio source and Be/X-ray binary LSI 61° 303 (Gre- 



gory & Taylor 1978; for recent work see Strickman 
et al. 1998). It has been suggested (Maraschi & 
Treves 1981) that the source contains a young pulsar 
and that the gamma-rays arise from inverse-Compton 
scattering of the optical photons from the B-star in 
the relativistic wind from the pulsar at the bound- 
ary between the pulsar wind and the stellar wind. 
However, 2CG 135-1-01 is a variable gamma-ray source 
(Tavani et al. 1998). If 2EG J0635-F0521 is similar 
to 2CG 135+01, then the constancy of its gamma-ray 
fiux (McLaughlin et al. 1996) would suggest a low 
eccentricity orbit. 

Finally, it may be that the emission is powered 
not by rotation of a neutron star, but rather by 
accretion. The recently reported detection of GeV 
pulsations from Cen X-3 with a frequency equal to 
that detected contemporaneously in x-rays (Vestrand, 
Sreekumar, & Mori 1997), albeit not with an out- 
standingly high significance, suggests the possibility 
that 2EG J0635-I-0521 could be an accretion powered 
system. The mechanism for gamma-ray production 
could be shock acceleration or a magnetospheric pro- 
cess. If SAX J0635-t-0533 is similar to Cen X-3, the 
x-ray emission would be from accretion and should 
be pulsed. However, the lower bound we obtained for 
the exponential cutoff energy of the x-ray spectrum 
is significantly higher than the cutoffs typically found 
in high-mass x-ray binaries (White, Swank, & Holt 
1983) and would be more suggestive of either a low 
magnetic field neutron star or a black hole. 

Identification of SAX J0635-I-0533 with the gamma- 
ray source 2EG J0635-I-0521 would argue against in- 
terpreting the gamma-ray emission as due to cosmic- 
ray production via shock acceleration in the super- 
nova remnant (Esposito et al. 1996; Jaffe et al. 1997). 
We note that we detect no diffuse hard x-ray emission 
in the BeppoSAX data. 

Possible ways to strength the identification of the 
x-ray /optical source and the gamma-ray source would 
be to detect a periodicity in the x-ray, optical, or radio 
which could then be detected in gamma-rays or to ex- 
tend the x-ray spectrum and join it to the gamma-ray 
spectrum. Finally, if the gamma-ray emission extends 
to higher energies, an improved gamma-ray position 
may be obtained from a ground-based TeV telescope 
such as the Whipple. 

We thank Fabrizio Fiore and the Sax Data Center 
for assistance with the SAX data. We thank Arlin 
Crotts for obtaining identification spectra at MDM 
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Table 1 

Lines Measured from the Spectrum 



Feature 


EW 


Flux (erg 


:cm "^s ^) 




FWHM 




(A) 


Observed 


Dereddened'^ 


(kms-1) 




Emission Features 








H/3 


3.22 


6.69 X 10-'^'^ 


3.91 X 10" 


-12 


400 


He I A5876 


1.31 


4.67 X 10"" 


1.13 X 10" 


-12 


600 


Ha 


30.3 


8.89 X 10-13 


1.39 X 10" 


-11 


390 


He I A6678 


0.40 


1.18 X 10-^^ 


1.74 X 10" 


-13 


450 


He I A7064 


0.84 


2.57 X 10"" 


3.17 X 10- 


-13 


470 


I A8446, Pais 


3.35 


8.94 X 10"" 


6.39 X 10- 


-13 


490 


Pal7 


2.10 


5.51 X 10"" 


3.89 X 10- 


-13 


570 


Pal6 


2.72 


7.08 X 10"" 


4.88 X 10- 


-13 


490 


Pal5 


1.73 


4.51 X 10"" 


3.03 X 10- 


-13 


470 




Absorption 


Features 








IS. A4430 


2.9 ±0.1 










I.S. AA5778,5780,5797 2.5 ± 0.2 










Na I AA5890,5896 


2.1 ±0.2 










I.S. A6283 


2.0 ±0.1 











Assuming E{B - V) = 1.2. 
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